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I. Introduction 



Interest in non-Abelian Einstein- Yang-Mills theories was sparked by the 
discovery of particle- likJ and non-Reissner- Nordstrom ( "hairy" ) black holeB 
solutions when the gauge group is su(2). Since then a plethora of hairy 
black holes, possessing non-trivial geometry and field structure outside the 
event horizon, have been found (see, for example,!), including coloured black 
holes in su{N) Einstein- Yang-Mills theories.&i Many of these objects have 
as a fundamental requirement for the existence of hair, a non-Abelian gauge 
field, often coupled to other fields, such as a Higgs field.! As in,i the existence 
of gauge hair is not surprising in itself, since the gauge field force is long- 
range. However, the non-Abelian nature of the field is important for evading 
the no-hair theorem, for example, for scalar fields coupled to the gauge field. 
The vast majority of these solutions have been found only numerically, with 
analytic work in this area at present being limited to an extensive study of 
the su(2) case,@'iil and some analysis of the field equations for general 

In this paper we continue this analysis of the coupled Einstein- Yang-Mills 
equations for an su{N) gauge field, and prove analytically the existence of 
"genuine" hairy black hole solutions for every A^. By a "genuine" su{N) 
black hole, we mean a solution which is not simply the result of embedding a 
smaller gauge group in su(A^). As in the su(2) case, the solutions are labeled 



by the number of nodes of each of the — 1 non-zero functions required 
to describe the gauge field. We shall prove that for each integer ni, there 
are an infinite number of sequences of integers rijv-i > ^^-2 > • • • > "^i 
corresponding to black hole solutions. It is to be expected that in fact every 
such sequence corresponds to a black hole solution, but unfortunately we are 
unable to prove this analytically, although we shall present a numerically- 
based argument for 5u(3) and numerical investigations (such as that done for 
su(5) in@) for higher dimension groups which indicate that this is in fact the 
case. The method used to prove the main theorem of this paper is remarkably 
simple, drawing only on elementary topological ideas. 

The structure of the paper is as follows. In section 2 we review briefly the 
su(A^) Einstein- Yang-Mills field equations and the ansatz and notations we 
shall employ in the rest of the paper. Next we state some elementary proper- 
ties of these equations, including results from.0 The remainder of the paper 
has a similar progression of ideas as,0 and we continue by first analyzing the 
behaviour of solutions to the field equations in the two asymptotic regimes, at 
infinity and close to the event horizon. Although these forms were discussed 
in,0 we present here a shorter proof. Section 5 is devoted to a discussion of 
the flat space solutions, which will be important for later propositions. The 
results here are somewhat weaker than in@ due to the fact that for su(A^), we 
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have — 1 variables and therefore have a 2(A^ — l)-dimensional phase space 
rather than a phase plane as in the su(2) case, so the powerful Poincare- 
Bendixson theory no longer applies. We now consider integrating the field 
equations outward from the event horizon and consider the various possible 
behaviours of the resulting solutions. As in,i there are three types of solu- 
tion: the regular black holes we are seeking, singular solutions (in which the 
lapse function vanishes outside the event horizon), and oscillating solutions 
in which the geometry is not asymptotically flat. The main results of this 
paper are in sections 7 and 8. An inductive argument is used to prove the 
existence of solutions for su{N) assuming existence for su(A^ — 1) (since we 
have rigorous theorems for su(2)i'i). The argument is presented in detail for 
su(3) in section 7 and a brief outline of the extension to general A^ in section 
8. Finally, a summary and our conclusions are presented in section 9. 

II. Ansatz and field equations 

In this section we first describe the ansatz we are using and outline the 
field equations. 

The field equations for an 5u{N) Yang- Mills gauge field coupled to gravity 
have been derived i n0 for a spherically symmetric geometry. We take the 
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line element, in the usual Schwarzschild co-ordinates, to be 

ds^ = -S^^i de + /^-^ dr^ + de^ + sin^ ^ (1) 

where the metric functions /i and S are functions of r alone and 

, , 2m(r) , , 

/x(r) = 1 ^. 2 

r 

A spherically symmetric su(A^) gauge potential may be written in the form0 

A = Adt \ B dr \ ]^{C - C^) \ [(C + C^)sin^ + Z}cos^] # (3) 

where D = Diagj/ci, k2, ■ ■ ■ , /^at} with ki > k2 > ■ ■ ■ > kj\f integers whose 
sum is zero. In addition, C is a strictly upper triangular complex matrix 
such that Cij ^ only if fc, = kj + 2 and its Hermitian conjugate, and 
A, B are anti-Hermitian matrices that commute with D. An irreducible 
representation of su{N) can be constructed by taking0 

D = Diag {iV - 1, iV - 3, . . . , -A^ + 3, -A^ + 1} . (4) 

In this case A, B have trace zero and can be written as 



for real functions ^fc, and similarly for 5. The only non- vanishing entries of 
C are 

C,, ,+1 = u;,-e^^^ (6) 



where uj and are real functions. All the functions in this ansatz depend 
only on r. Now we make the following simplifying assumption^: 



(7) 



The remaining gauge freedom can then be used to set B = 0,113 which implies 
that the are constants, which we choose to be zero for simplicity. 



The gauge field equations for the tUj then take the formti!Hiil(where we 



have set K = 2 in0) 



1 - + 2 «i + 



0, (8) 



where 



Pe 



1 ^ 



(9) 



and the Einstein equations can be simplified to 



/ ^ 2 ^ 2Gr 

m = nG + r Pe, — = — 

b r 



(10) 



where 



N-l 



'11^ 



i=i 



Note that ujq = = ujn, so that these are the usual equations! in the su(2) 
case, and have a very similar, but slightly coupled, structure for general A^. 
The equations (PHTTD possess two symmetries.0 Firstly, the substitution 



—ujj for any fixed j leaves the equations invariant, exactly as in the 



su(2) case. Secondly, there is an additional symmetry under the transforma- 
tion j ^ N — j for all j. We will not assume that this symmetry is respected 
by the solutions of the field equations, so that the 5u(3) case is not necessarily 
trivial. 

We are concerned in this paper with black holes, which will have a regular 
event horizon at r = rh, where /i = 0. The equations in their present form 
are singular at r^, so in order to produce a set of equations which are regular 
as /i ^ 0, we define a new independent variable r by0 

dr ^ /X 
Tr = <^2) 

denote d/dr by', and define new dependent variables k, Uj and \1/ as followfB: 

^ = ^, U, = ^uj'^, «: = ^ (1 + + - 2r2p,) . (13) 

Then the field equations take the form!: 

r = r\l/ (14) 

Uj = rUj (15) 

^ = (K-^)^-2/iG (16) 

(5^)- = 5(k-^)^ (17) 

LOj (18) 

k = + l + 2fxG. (19) 



The main thrust of this article is to prove the existence of regular, black 
hole, solutions of the field equations for general with — 1 degrees of 
freedom. In other words, we begin with a regular event horizon at r = r/j, 
where /i = 0, and integrate outwards with increasing r. Later we shall classify 
the possible behaviour of the solutions as r increases. We are interested 
primarily in those solutions which possess the physical properties of a black 
hole geometry, namely for which fi > for all r > rh, the ujj and their 
derivatives are finite for all r > and the spacetime becomes flat in the 
limit r oo. We shall refer to such solutions as regular black hole solutions. 

III. Elementary results 

In this this section we state a few elementary results, which will prove 
useful in later analysis. Firstly, two lemmas which are proved in the su(iV) 
case exactly as in the su(2) case, see.i 

Lemma III.l // /i(ro) < 1 for some vq then fi{r) < 1 for all r > r^. 

Lemma III. 2 As long as < fi < 1 all field variables are regular functions 
ofr. 

These two lemmas show that for a regular event horizon at r = rh, where 
= 0, then /i < 1 for all r > r^ and the field variables are regular functions as 
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long as /i > 0. The black hole solutions we seek approach the Schwarzschild 
geometry as r — > oo, so we may assume that all ujj are bounded for all r 



(from lemma [1II.2| each ujj is bounded on every closed interval [r/^, ri], so we 



are simply assuming that ujj remains bounded as r oo). It will be proved 



in section that this assumption is in fact valid for solutions having /i > 



for all r. Define a quantity M, for each j to be the lowest upper bound on 



ujI i.e. 



^|<M,, Vj = l,2,...,iV-l. (20) 

It has been shown ini that the following result is true for N = 2, and proved 
in general with the assumption that all uj are bounded via an elegant method 
in.0 

Theorem III.3 Mj < j{N - j). 

Part of the proof of this theorem in0 involves a knowledge of where ujj 
may have maxima or minima. The result following is less powerful than in 
the su(2) case due to the coupling in the gauge field equation (||). 

Proposition III. 4 As long as fi{r) > 0, the function ujj{r) cannot have 
maxima in the regions 



1 . o n X . „ L 1 



(21) 



or minima in the regions 



and < < ^/ 1 + - + . 

(22) 



Proof 



When = 0, from 



J' 



(23) 



so that ujj will have a maximum if 



ujj > and u;| < 1 + - (t^j+i + 



(24) 



or 



ujj < and > 1 + - + ^^j-i) 



(25) 



Similarly, Uj will have a minimum if 



> and ujj>l + - (t^j+i + ^"j-i) 



(26) 



or 



< and < 1 + - + ^"j-i) 



(27) 
□ 
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IV. Asymptotic behaviour of solutions 

We seek solutions to the above field equations representing black holes 
with a regular event horizon at r = and finite total energy density. In 
order to prove the local existence of solutions of the field equations with the 
desired asymptotic behaviour, we shall apply the following theorem.0 

Theorem IV. 1 Consider a system of differential equations for n + m func- 
tions a = {ai,a2, ■ ■ . , an) and b = {bi, 62, • • • , &m) of the form 

= xP'fi{x,a,b), 

ax 
db- 

x-^ = -Xibi + x''^g,{x,a,b), (28) 
ax 

with constants Aj > and integers Pi, Qi > 1 and let C be an open subset 
of such that the functions fi and Qi are analytic in a neighbourhood of 
X = 0, a = c, b = 0, for all c E C. Then there exists an n-parameter family 



of solutions of the system (jl^ such that 



a.,{x) = Q + O(x^0, H^) = (29) 

where ai{x) and bi{x) are defined for c E C, and \x\ < Xq{c) and are analytic 
in X and c. 



In this section we consider only the equations for /i and the ujj. The equa- 
tion for S (^) is independent of the others and hence can be integrated 
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immediately once we have the other field variables, with an additional pa- 
rameter. This procedure gives a finite answer on any closed, bounded interval 
on which the remaining variables are finite. Results similar to propositions 
IV. 2| and |IV.3| were proved in,0 but here we are using a more straightforward 



approach. 



A. Behaviour at infinity 

As r — > oo, in order to have finite total energy density, the geometry must 
approach the Schwarzschild solution, that is 

m — ^ M = constant, G ^ 0, pg 0. (30) 

Hence each ujj must approach a constant value, which is fixed by (^Up to be 

co^^jiN-j). (31) 

The equations for the Uj are then automatically satisfied in the limit as 
r — >• oo. 

Proposition IV. 2 There exists an N -parameter family of local solutions of 
(]^-[77y near r = oo analytic in Cj, M and such that 

^.(r) = VjW^)-'f + 0(^-^y (32) 
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Proof 

Introduce new variables 

= ipj (^y/j{N - j) -ujj^ , ^j^r'^iiuj'j A = r(l-//). 

(33) 



X 



Then 



where 



N-l ^2 

G = -'E| (35) 



and 



AT 



that is, 



+2xV',-iV(j-l)(^-J + l)l ' , (36) 



where is a polynomial in x, , ^'s and ■^'s. Similarly, 

where j^. is analytic in x. A, ^'s and ^'s in a neighbourhood of x = 0. Using 
the field equations we also have 

= -2^,- + - j)] - V^.-i - 1)(A^ - j)(^ - J + 1) 

-^j+iVJU + 1){N- j){N - J - 1) + a;4 (39) 
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where f^. is a polynomial in x, n ^, ■^'s and ^'s. The algebra simplifies if we 
define new functions aj, (5j such that 

= a,- + A/3,- (40) 

where A is a real constant not equal to unity. The equations for and (5j 
then read 

(A-l)x^ = 2a,- - 2\j{N - j)]a,- + Vi(j - 1)(A^ - j)(^ - J + 

+ + (A' + A)/?,- 

-2b-(iV - j)]/?i + ^j{j-l){N-3){N-3 + l)P^_, 

(A-l)x^ = -2a,- + 2[7-(A^ - j)]«,- - - 1)(A^ - 3){N - j + l)a,_i 

-,/j{jTW^^jW^r^<^J+i + (1 - 3A)/3,- 
+2b-(iV - - ^j(j-l){N-j){N-j + l)Pj_, 
-Vjij + ^)iN-j){N-j-l)P,+, + xff,^, (41) 
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where the fa^ and //j^ are analytic functions of x, ^, a's and /5's. Consider 
the matrix A^^v-i whose entries are 

/ „ ™ . \ 



2(A^-1) -^(A^- l)2(iV- 2) 

-^y{N -1)2{N -2) 2.2(iV-2) 



-^2{N -2)3{N ~3) 



J 



(42) 



For a vector q = {q^, . . . , qn-i), let pj = \/ 3{N - j)qj for j = 1, . . . , - 1. 
Then 



q MN-iq = Pi+PN-i + iPi-P2) + . . .{pn-2- Pn-i) ■ (43) 

Hence the matrix Ai^^i is real, symmetric and positive definite, and will 
have positive real eigenvalues Si, and corresponding eigenvectors Vi. We may 
expand our variable vectors in terms of eigenvectors of Ai^-i as follows: 

N~l N-l 

i=l i=l 

N-l N-l 



(45) 



=1 



i=l 



in terms of which the equations (^T]) now become: 



(A-l)x 
(A - l)x 



dAi 
dx 
dB^ 
dx 



(2 - £i)Ai + (A^ + A - £i)Bi + xFi 
{£i - 2)Ai + (1 - 3A + £i)Bi + xd. 



(46) 
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This transformation has removed the couphng between the j's but the A's 
and -B's are still coupled. Repeating the procedure with the matrix Vi will 



decouple these equations, where 

/ 



\ 



2-S, + A - 
-2 + £i 1 - 3A + 



This matrix has eigenvalues 



/x± = ^(l-A) (3±v^r+4^ 



and corresponding eigenvectors V±. Writing 



the equations are now 



(47) 



(48) 



D_(x)V_, 
(49) 



dC 

(A - l)x—^ = fi±C± + xD±. (50) 
dx 



In the case where < 0, theorem |1V.1| applies directly and we have 



C_ = 0(x); (51) 



when /i_ = (corresponding to £i = 2), 



C^ = K_ + 0{x) 



(52) 



for some constant K_. The other eigenvalue //+ is always strictly positive, 
and in this situation theorem [1V.1| does not apply. However, equation (|50| ) 
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can be integrated using the standard Picard method (see, for example, to 
give an analytic solution 

C+ = K+ + 0{x). (53) 

Both C± will be analytic in K± and x in some neighbourhood of a; = and 
substituting back through equations (0,^) gives solutions of the form (|32| ) 
as required. □ 

B. Behaviour at the event horizon 

We assume that aXr = Vh there is a non-degenerate event horizon, namely 
Ai(r/i) = 0, but n'irh) > is finite. 

Proposition IV. 3 There exists an N -parameter family of local solutions of 
(^-[77|j near r = r^, analytic in r^, Uj^h and r such that 

K^h + p) = fi'{rh) + 0{p), 

^j{rh + p) = ujj,h + oo'jirh) + 0{p^) (54) 
where p'{rh) and oj'jirh) are functions of the Uj^h- 
Proof 

Let p = r — rh he the new independent variable, and define 



X = r, A = — , i^j = ^j, = ^— ^- (55) 
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Then the field equations take the form 

dx 
^ dp 
dX 



P 



P 



dp 

dipj 
dp 
d^ 
dp 



-A + 



(56) 



1 



X 



- "^xpe 



+ P- [1-2G] 

X 



= -X + pHx + Fx 

P§i 
X 

P§^ = -^J + PH,+F„ 



(57) 
(58) 
(59) 



where the F's and H's are polynomials in a; ^, A ^, and the other variables, 
and the F's depend only on x and ipj^s. Next define 



Cj-Cj-^j^ X-X-Fx 



(60) 



whose derivatives are given by 

dp 
dA 

dp 



P 
P 



<j+pGj 
-X + pGx. 



(61) 
(62) 



Here the G"s are analytic in x ^, A ^, X, x, A, ipj, 0- ^PPlyi^g theorem [IV. 1 
there exist solutions of the form 



X 



rh + P, ipj = ^j,h + O (p) , X,ij = (p) , 



(63) 



which gives the behaviour (0), together with the required analyticity. From 
the field equations (|^-[TT|), setting p{rh) = gives the following relations: 

1 



P\rh) = 7:rhPeirh) 



(64) 
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where 



N 



Pe[rh) 



N 



2iT 



(66) 
□ 



V. Flat space solutions 

The behaviour of the gauge field equations in the fiat space hmit will be 



useful later in section ^ when we come to examine the properties of regular 
solutions. We are not concerned here with the global existence of fiat space 
solutions but only the local properties pertinent to the curved space problem. 

The analysis here is similar to that in,i but note that for general N there 
will be A^— 1 coupled gauge field equations. The powerful Poincare-Bendixson 
theory of autonomous systems is not applicable when N > 2, and so we will 
be able to derive only correspondingly weaker results, and cannot draw a 
phase portrait. Fortunately the theorems later require only a knowledge of 
the nature of the critical points and the local behaviour close to these points, 
which is the subject of this section. 

In fiat space the gauge field equations reduce to: 

= 0. (67) 
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2 // I 

r uj^ + 



These equations can be made autonomous by changing variables to r = log r, 
and denoting d/dr by ': 



ib^ — UJj + 



uj-i — 0. 



This system of — 1 coupled equations has critical points when 



l-^l+ok'+l+^l-l) 



If LOj 7^ for all j then 



cu.^O. Vj = l,2,...,7V-l. 



= Vj{N-j). 



There is also a critical point when cuj = for all j. 



(68) 



(69) 



(70) 



The other critical points can be described as follows. Suppose that cui — 
and uji+k = but that uJi+m 7^ for m = 1, . . . , k — 1 (where k > 2). The 
critical point is then described by the equations 







= 1 



2 I / 2 I 2 \ 



— u. 



+k 



(71) 



which have the solution 



\/ m{k — m) lor m— 1, . . . , /c — 1. 



(72) 
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In other words, if we have a run of k non-zero a;'s with zero Ui at each end, 
then the solution for the non-zero a;'s is the same as the N — k case with no 
zero a;'s. We can also put together a run of as many zero cu's as we like in 
the critical point. 

In order to classify the critical points, we linearize the field equations. 

Let 



where ouj is the value of ouj at the critical point and Sj is a small perturbation. 
The equation for ej is, to first order. 




(73) 








_L, ,(0) 0. , ,(0) _L . , ,(0) _L . , ,(0) 



(74) 



We shall take the two cases we need to consider in turn. 



Firstly, the case where a;] — 0, when the equation for reduces to 




(75) 



In order to find the nature of the critical point, let 




(76) 



then A satisfies the equation 




(77) 
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This implies that 



1 1 



where a is the positive real number given by 

a' = 3 + 2u;fll + 2u;l^l (79) 

We conclude that, in the plane, we have an unstable focal point, as 

found in the su(2) case in.0 



Secondly, we consider the situation in which there is at least one ci;j°^ 
which is non-zero. Suppose that = 0, and ool'^j^ = 0, but ool'^j^ ^ for 
m = 1, . . . , k — 1, where we include the case that i = and k = N. Then we 
have a series of coupled perturbation equations: 



n _ •■ • o (0)2 I (0) (0) 



n _ •■ • o (0)2 I (0) (0) /Qf.N 

U — Cj+fe-l — Q+fc-l — ^ei+fc-l^^j+jfc-i + ei+A:-2^^j+fc-2^i+A:-l- ) 

Define a vector e by 

e = (cj+i, • • • ! ej+yfc-i)^ (81) 
and consider solutions of the form 



e = e^^qr (82) 
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where g is a constant vector. Then A^ — A are eigenvalues of the matrix Aik-i 
given by (|^). As discussed in section the matrix Aik~i has positive real 
eigenvalues £j. Then 

A = \ ^=j + l,-j (83) 

will have positive and negative real values and there is a saddle point. Again 
this is in direct analogy with the su(2) case.i 



VI. Global behaviour of the solutions 

In this section we investigate the behaviour of solutions of the regular 
field equations (|T^-|I^) as functions of r. From section we know that, 
given any starting values iOj^h for the gauge field functions, then there is a 
local solution of the field equations in a neighbourhood of the regular event 
horizon at r = r/^, r = 0, which is analytic in r and the initial parameters. 



Furthermore, from lemma |III.2| as long as < \E' < 1 the solutions remain 
regular. Therefore, as we integrate out in r from the event horizon there are 
only three possibilities: 

1. There is a tq > such that \E'(to) = 0. 

2. For all r > 0, we have \E'(r) > and r(r) remains bounded as r ^ oo. 

3. For all r > 0, we have ^'(t) > and r(r) ^ oo as r ^ oo. 
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We refer to solutions of the first type as singular solutions. Those of type 
2 are known as oscillating solutions ini and we retain their terminology. It 



will be shown below (proposition |VI.4|) , that there is a maximum value r 



of r for oscillating solutions. Finally, we denote by 5*00 solutions of the first 
type for which r(ro) < r^ax- Our first task is to show that solutions of type 
3 are precisely the black hole solutions we are seeking. We begin by proving 
the assumption made at the end of section ||, namely that for solutions of 
type 3, each uj is bounded for all r. 

Lemma VI. 1 // \1/(t) > for all t and limT-^oo ''"(t) = oo then all the Uj 
remain bounded as r ^ oo. 

Proof 



From the local existence propositions |IV.3| and |IV.2| and lemma |IL2|, each 



ujj is an analytic function of r as long as /x > 0. Introduce a new variable 
X = r~^, then each Uj is an analytic function of x in a neighbourhood of 
X = 0, except possibly at x = 0, and can be written as a Laurent series: 

oo 

u,{x) = J2 (84) 

n=— oo 

In order for > for all r, it must be the case that 2mr~^ = 2mx < 
1, although m itself need not necessarily remain bounded. Hence 2mx is 
bounded in a neighbourhood of x = 0, and in particular is analytic at x = 0. 
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Therefore we can write 



2mx = hnX^. (85) 



n=0 



Thus 



^ = 2 S ^-"^ ' ^ ""^ ^ " 2 S ^"^^ " ^^"^ ' ^^^^ 

n=0 n=0 

SO that dm/dr is analytic in a neighbourhood of x = 0. Prom the field 
equations, 

dm ^ , s. 

— = I^G + r^Pe (87) 

where 

N-l / , N 2 2 ^ 

G=E(^ r=p. = ^E[c.J-a,^,-iV-l + 2,f. (88) 

Since both G and r'^pe are positive, they must each be analytic in a neigh- 
bourhood of x = 0. Consider G first, using 

^ = -x^^ = - y ainx^^' (89) 
dr dx ^ ^ ' 

n=—oo 

and since must be analytic near x = 0, it must be the case that 

al = Vn < -1. (90) 

Then 



n=-l 
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for constants and hence 

a;| - a;|_, - iV - 1 + 2j = - [a^\ - a^^^'^'] + J] (92) 

n=— 1 

Since r^pe is analytic, 

ail = ±«r/ Vj. (93) 

But loq = 0, so aP_i — and thus a'_^ ~ for all j. We conclude that Uj is 
analj^ic at x = 0, and therefore finite as x = 0. Therefore ujj is bounded for 
all r e [r/j, oo). □ 

Proposition VI. 2 If ^{t) > /or all r and lim^_^oo '^('T') = ^^en ^/le 

solution tends to one of the fiat space critical points, with the exception of 
the origin. 

The proof of this proposition closely follows Proposition 14 of Ref. 7, the 
main thrust of which is the following lemma. 

Lemma VI. 3 //^'(t) > for all r and lim^^oo r(r) = oo then 

lim *(t) = 1. (94) 

T— ►oo 

Proof 

There are two cases to consider. 

1. There is at least one cuj which has zeros for arbitrarily large r. 
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2. All uij have only a finite number of zeros. 



Case 1 

The proof in this situation follows exactly that of Proposition 14 of Ref. 7 
(equations (74-76)) since all the cuj's are bounded (lemma |VI.1|) . 



Case 2 

In this case there is a Ti > such that for all r > Ti no ujj has a zero. 
However, since the equations governing the Uj are coupled, the ojj do not 



necessarily have to be monotonic, provided proposition III. 4 is satisfied 



Suppose, initially, that there is a T2 > Ti such that for all t > T2 every 
Uj is monotonic. In this situation each uj has a limit and so limr-,00 Uj = 
for all j. Then the proof of Proposition 14 of Ref. 7 carries over directly to 
show that m(r) is bounded and limT-_,oo ^(t) = 1. The proof proceeds as 
follows. Firstly, in analogy to equation (74) of Ref. 7, integrating the field 
equations (P^|-[T^) gives, for each j, 

|W,(r2)-W,(n)|<^ (95) 

for all T2 > Ti > T2, and some constant Cj, since all the lvj are bounded. Fix 
Ti for the moment, then for all > Ti, 

r rmf dr' < Cj r Uj dr' < c'- (96) 
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for constants Cj- and c^', since uj has a limit. Hence 



N-l 

m{r,)-m{n)<-Y,c; (97) 

which is bounded. 

Next turn to the other extreme situation, where has minima for arbi- 
trarily large r, for every j. If Uj has a minimum at r = tq, then by proposition 



^J(ro) > 1 + ^ (^J+i(ro) + <i(ro)) • (98) 

Define Nj to be the greatest lower bound of the set of minimum values of u'^ 
for r > Ti, where Nj can be zero (although uJj cannot). Then we have 

iV, >l + i(iV,+i + iV,_i). (99) 

These inequalities may be solved to give 

N,>j{N-j) Wj. (100) 



However, since Uj < j{N — j) for all j from theorem |III.3| , it follows that 
uJj{T) = j{N — j) for sufficiently large r and all j. Therefore Uj{T) is zero for 
sufficiently large r, hence m(r) is bounded as r — >■ oo and limT-_^oo ^{'t) = 1- 

The remaining intermediate case is where ujj and cUj+j, where / > 1, are 
monotonic for all sufficiently large r, whilst z = 1, — 1 have minima 
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for arbitrarily large r. We include the possibility that j = 0, to cover the 
case where there is only one Ui which has a limit. Define Li and L2 by 

uj] Li, uj]_,i ^ L2 (101) 

as r — s> 00. Given e > 0, there is a Ti such that 



\uj]-Li\<e, \uj]^i-L2\<e Vr > Ti. (102) 



If ujj_^i has a minimum at Tq > Ti, then 



o^|+i(ro) > 1 + ^ (^1+2(^0) +^|(ro)) 



2 
1 
2 

With A^,+i as before, then 



> 1 + 1 (^2^2(^0) + ^! -e). (103) 



iV,+i > 1 + + (104) 

where Ki = ^ (Li — e). Similarly, 

iV,+/-i > 1 + + K2 (105) 

where K2 = \ {L2 — e). The remaining inequalities for Nj^2, ■ ■ ■ , ^j+1-2 are 
exactly as before (P5|). The previous method can be repeated to give 



Nj+i >i{l-i) + Ki + K2 z = 1, ...,/- 1. (106) 

Now define iVj to be the lowest upper bound of the set of maximum values 
of (which exists since ujf is bounded). If uJj_^_'^ has a maximum at tq > Ti, 
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then 

a;|+i(ro) > l + i(a;|+2(To)+a;|(ro)) 

> l + ^{u;]^, + L, + e). (107) 

The analysis now proceeds exactly as in the previous paragraph and yields 

Nj+i < i{l - i) + Ki + K2 (108) 

where Ki — ^ {Li + e) for i = 1, 2. In addition, by definition it must be the 
case that 

Nj+i > N^M (109) 

which means that 

i{l - ^) + i (Li + L2) - e < Nj+i < Nj+i < i(/ - ^) + 1 (Li + L2) + e (110) 
for all e > 0, whence Nj+i = Nj+i for all i and 

u^]+i^i{l-i) + \{L, + L2) (111) 

for sufficiently large r. 

In conclusion, then, we have shown that the intermediate case has some 
a;j's which are monotonic for sufficiently large r, with the remaining ujj''s 
being constant for sufficiently large r. Hence in this case also m is bounded 
as T — > 00 and lim^^oo ^{^) = 1- D 
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Proof of Proposition |V1.2| 



In order to show that the geometry becomes flat as r — > oo, it remains to 
show that 5 ^ 1 as r ^ oo. From the proof of lemma |V1.3| , for ujj having a 
limit as r ^ oo, we have 

rUi = dji^O (112) 



as r ^ oo. In case 1 of the proof of lemma |VI.3| , the proof ol0 carries straight 



over to show that ?7j ^ as r ^ oo also in this situation. Now 

N-l 



1=1 



for some constant C and sufficiently large r. Therefore 

_ SJ_ ^ C ^ C 
S S r\E' ~ fir^ ~ 

for some constant C and sufficiently large t, as fi 1 as t 
means that S has a finite limit as r — oo. 



:ii3) 



(114) 
oo. This 



The field equations only involve ^ rather than just 5* itself. Therefore S 
is defined only up to a multiplicative constant, and without loss of generality 
we may therefore take the finite limit of S" as r ^ oo to be 1, so that the 
spacetime is asymptotically fiat. 

Let 6{t) = 2\E' — k — 1 be a small perturbation, then the equation for ujj 



reads 



Uj + 



UJj = {l + 5{T))Uj 



(115) 
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as r ^ oo. Since 6 is small, it does not alter the position or nature of the 
critical points as compared with the exactly flat space case. Lemma |Vl.ci 



showed that each djj,ujj as r — > oo. Hence Uj must approach one of 
the critical points whose nature was elucidated in section 0. The flat space 
analysis showed that, if Uj = 0, then there is an unstable focal point in the 
{uj, ujj) plane. Hence our solution cannot approach this value, unless ujj = 0. 
The solution has to tend to one of the saddle points along a stable direction, 
and, in direct analogy with the su(2) case,lll there are no solutions where 
— > as r — s> oo. The solution must therefore be a member of the family 
found in the local existence theorem |1V.2|. □ 



The power of proposition [VI. 2| lies in that, if we can prove the existence 
of solutions for which \E' > for all r > and r — > oo, then these solutions 
are automatically the regular black holes (and, if r/^ — 0, soliton solutions^) 
we seek. We close this section by determining the asymptotic behaviour of 
solutions of type 2. 

Proposition VI. 4 // \E'(r) > for all t > and r(r) remains bounded, 
then ^ and n —>■ 1 as t oo. In addition, there is at least one j for 
which ujj as T —>■ oo, and this ujj has infinitely many zeros. 

Proof 
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Since r is monotonic increasing and bounded, it has a limit tq. In addition, 
m is monotonic increasing and bounded since the positivity of fi imphes that 



m < - < — . 
2-2 



(116) 



As r — s> oo, r ^ and hence \1/ — > from ([T^)- Consider the quantity 

2 2 

= - J (1 + - 2/t^) = (-^G + r^pe) • (117) 



Then 



< r^/iG'(l - k) 



< 0, 



:ii8) 



for sufficiently large r for which \E' < 1/2, and since k > l.i Hence 



4 



:ii9) 



as r ^ oo, since E is monotonically decreasing for sufficiently large r. Then, 
following,0 we have k — 1 as r — * oo. 

Let 6{t) = k, — 2\E' — 1 be small, then the equation for uj reads 



UJj = UJj 



(120) 



Replacing r by — r yields the equation 



:i2r 
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which is the equation for flat space solutions, so that uoj must approach one 
of the critical points. The fact that 7^ means that pq cannot tend to zero 
as r — > 00 because ixG does vanish in the limit. Hence at least one of the 
cjj's must be zero as r — 00. Since 5 is small it alters neither the position 
nor the characteristics of the critical points, and hence this uoj will go into 
the focus at zero. Therefore it oscillates infinitely many times before it hits 
zero. □ 

We may determine the value of tq as follows. Since k is finite and we have 
( |119| ), it follows that r'^pe — s>l/2asr^oo. Hence 



The maximum possible value of Tq is when all cjj ^ as r ^ 00. In this 
case: 




(122) 




(123) 



The minimum possible value of Tq is when only one ujj ^ as r 



cx), when 



(124) 



which has a minimum when j = 1 or — 1, and 




{N-lf. 



(125) 



With this value of rg, we 



and denote by S, 



00 



singular solutions for 



which jj, vanishes outside the event horizon when r < Tq, 



.max ■ 
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VII. su(3) black holes 

We are now in a position to prove the existence of regular black holes. We 
begin in this section with su(3), which we shall study in some detail, before 
proceeding to the general case in the next section. The method in general 
is exactly the same as in this section, but it is hoped that by doing the 
simplest specific case, where we have just two gauge field degrees of freedom, 
explicitly, the proof will be more transparent. 

Genuine su(3) solutions (as opposed to embedded 5u(2) solutions which 
are discussed below) have been found numerically in.0 There it was con- 
jectured that, analogous to the su(2) case, black hole solutions exist for an 
infinite, but discrete, set of points in the parameter space given by the values 
at the event horizon of the functions describing the gauge degrees of freedom. 

Firstly, we review briefly the way in which su(2) solutions may be embed- 
ded in su(3). For neutral black holes, numerical solutions are discussed in,i 
where both su(2) and so (3) embeddings are discussed. The so (3) embedding 
contains scaled su(2) black holes as well as genuine so (3) solutions possessing 
two degrees of freedom, which are indexed by the number of nodes of each 
of the gauge field functions. The method for embedding charged black holes 
can be found in,l where the various possible embeddings are described in de- 
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tail and numerical solutions presented for the various embeddings in su(5). 
Such solutions possess Coulomb charge arising from one or more of the gauge 
degrees of freedom. The remaining gauge field functions do not contribute 
to the charge and have structures similar to the neutral solutions. The black 
holes are once again indexed by the number of nodes of the gauge functions. 

Neutral black holes 

Let uji{r) = v^co'(r) and t02ir) = \/2uj{r). Then the field equations ([^-pJ]) 
become: 

r^fxu" = -{2m-2r%)uj' - {l-uj^)uj, (126) 

m' = {fiG + r^pe), (127) 

S' _ 2G 
'S ~ 
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where 



/x = l-^, G = Au'^ pe = ^^{co'-l)\ (129) 

In order to extract the precise form of the su(2) equations, define a new 
independent variable R and a new dependent variable M by 

R = ^r, M = ^m, (130) 

with the other field variables remaining unchanged. Then the equations are 



2 (PuJ 



2M 



2R 
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^-(l-.V. (131) 



dM 
ItR 

ldS_ 
SdR 



f duo\ 1 . 9 



2 fdoj 
R [dR 



(132) 
(133) 



These are now exactly the su(2) field equations, and the existence of regular 
black hole solutions has been proved inffl 

Charged black holes 

Suppose uj2{r) = (the case c<Ji(r) = is identical). Then 



m 



(2m - 2r^po) uj[ - [l - Luj) lui, 



2J^ 



(134) 
(135) 
(136) 



While these equations are not the same as those for su(2) neutral black holes, 
the equations for k, \E', U, uji, r and S as functions of r are identical. Hence 
the analysis ofi carries over immediately to prove the existence of black hole 
solutions of this form. 

Our first existence proposition states the existence of genuine su(3) black 
holes, that is, black holes which do not belong to one of the two families 
outlined above. The proof is remarkably simple once we have set up our 



notation suitably. 



37 



Proposition VII. 1 There exist regular black hole solutions of the su(3) 
Einstein- Yang- Mills equations which do not correspond to embedded su(2) 
charged or neutral black holes. 

The strategy in proving this proposition is very similar to that used in.0 
Generic values of the starting parameters {uJi^h,^2,h) lead to a singular so- 
lution, but we will prove that there must be some values of the starting 
parameters which do not lead to a singular solution. Firstly, we rule out the 



possibility of solutions of type 2. From section any solution for which 
/i(r) > for all r and r(r) remains bounded as r oo, has r — >■ ro as 
r — > oo, where 

(A^-l)2 <r^ < 1a^(A^ + 1)(A^-1). (137) 

For = 3, therefore, tq = 2. Fixing > 2 will therefore rule out the 
possibility of such solutions. This is analogous to the value > 1 which 
rules out such solutions in the su(2) case. At the end of this section we shall 
return to the existence proof for <2. 

For every pair of starting values {uji^h,^^2,h) there are then just two pos- 
sibilities: 

1. /u(t) > for all r and we have a regular black hole solution; 

2. there is a tq > such that /x(ro) = and we have a singular solution. 
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Define a new variable R by 



R='-^ (138) 



and Rm by the maximum value of R for each solution, that is, in case one 
Rm = 1 (corresponding to r — >• oo), and in case two Rm = -R(to)- For each 
solution we define to be the number of zeros of the function Ui between 
r = and r = oo in case 1 and r = tq in case 2, where Ui can be infinite. 
Defining new variables Ni by 

AT, = (139) 
we allow the possibility that A'^j = 1. 

Each pair of non-zero starting values {uji^h,uj2,h) can then be mapped to 
the three quantities [Rm, Ni, N2) for the corresponding solution, giving rise 
to a map from ((f^i./i, ^^2,h) space) to x [0, 1] {{Rm, Ni, N2) space), where 
B is the discrete set {0, 1/2, 2/3, 3/4, . . . , 1}. Call this map /. Figures 1-3 
sketch the nature of {uJi^h,^2,h) space, B x [0,1] and B^ respectively. Note 
that this map will not be 1-1 for = 3, although it was conjectured i n@ that 
for N = 2 the map is 1-1. 

Note that since we know the nature of the solution space when one of the 
uji^h vanishes (since in this case Ui = 0), we do not need to extend the map / 
to the co-ordinate axes in {ui^h, ^2,h) space. With this notation, we are ready 

39 



to prove proposition |V11.1| , once we have the following lemma. 

Lemma VII. 2 Suppose that for starting values {uii^h, ^2,h) there is a singular 
solution with ^{tq) = and the gauge field functions having node structure 
(ni,?7,2)- Then the map f is continuous at {u!i,h,^2.h) ■ 

Proof 



From proposition |IV.3| , the field variables are continuous in r and the starting 
parameters. Thus all starting values {uji^h,^2,h) in a sufficiently small neigh- 
bourhood of {LJi^h,^2,h) will give rise to a singular solution with /i(ro) = 
where Tq is close to tq, the values of the field variables at tq will be close to 
those at tq in the original solution and the node structure will be (^1,^2) 



since the gauge field functions cannot have double zeros (proposition p.ll.4[) 



In other words, / is continuous at {uji^h, ^2,h) ■ D 



Proof of Proposition VII. 1 



Consider the open subset of the {uji^h,^2,h) plane given by 

D = {(tui,/,, uj2,h) ■■ < uji,h < uj2,h}- (140) 

The subset D' = {{uji^h,^^2,h) '■ < uj2,h < ^i,h} '^^^ t)e treated similarly. The 
symmetries of the field equations (pj-pJ]) mean that it is sufficient to consider 
only the first quadrant of the {uJi^h,^2,h) plane. Fromi we know that along 
the line uji^h = ^2,h there are singular solutions with node structure (ni,ni) 
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for all rii = 1,2,.... Therefore there are neighbourhoods of D corresponding 
to singular solutions with node structure (rii, ni) for all rii = 1, 2, . . .. Hence 
f{D) contains points corresponding to {Rm, Ni, Ni) for all Ni G B. Therefore 
f{D) is not a connected set. Therefore / cannot be continuous everywhere 
on D because D is connected. Hence we conclude that there exists at least 
one {ui^h,i^2,h) G D corresponding to a black hole solution because / is 
continuous for all {uji^h,^2,h) corresponding to singular solutions. □ 

This approach allows us to see quite simply how the transversality prop- 
erty conjectured in§ arises in the su(2) case. Along the line of values of Uh, 
singular solutions having different numbers of nodes must be separated by at 
least one regular solution. For su(2), it is known that values of ooh sufficiently 
close to 1 correspond to singular solutions with uj having one node, whilst 
sufficiently small Uh correspond to solutions having as many nodes as we 
like. In addition, if we have a regular solution with n nodes, then all singular 
solutions with Uh sufficiently close have either n or n + 1 nodes. Starting 
with the singular solutions with one node, and decreasing Uh, we first hit a 
regular solution with one node (there being only the trivial solution with no 
nodes). Then there may be more singular solutions with one node, or two 
nodes. In the former case there must be another regular solution with one 
node, in the latter case there will next be a regular solution with two nodes. 
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Either way, there must be a regular solution with two nodes before we can 
move on to singular solutions with n > 2 nodes. The process continues for 
all n. 

Having shown that there exist genuine su(3) black holes, the main result 
of this section is the following theorem, which we prove at the moment for 
the case rh > 2, returning to the case r/^ < 2 at the end of the section. 

Theorem VII. 3 Given fii = 0, 1, . . then there exist regular black hole 
solutions of the su(3) Einstein- Yang-Mills equations with uji{r) having fii 
nodes and uj2{r) having n2 nodes, for infinitely many n2 > ni- 

A similar result holds with the roles of ui and uj2 reversed. Note that this 
result is slightly weaker than the corresponding theorem for su(2), since we 
cannot guarantee that every combination of {ni,n2) is the node structure of 
the gauge fields for some black hole solution, only that an infinite number of 
such combinations does occur for each ni. At the end of this section we shall 
give an argument, based on a numerical analysis, that in fact black holes 
exist for all (77,1,^2), although we are not able to prove this analytically. The 
proof of theorem |VII.3| will proceed via a series of lemmas. 



Lemma VII. 4 Suppose that the starting parameters (0, uJ2,h) correspond to 
a charged regular black hole solution in which 0^2 1'^) has n2 zeros. Then, 
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given no, for all sufficiently small ooi^h ^nd uj2,h sufficiently close to 0O2,h, the 
solutions are regular or singular solutions with rii > tiq and 712 > ri2 . 

Proof 

By continuity, all field variables will remain close to the original charged 
solution until r ^ 1 and the geometry is approximately flat. At this point 
cui will be very small and CU2 will be close to 1. The equations for ei = cui 
and 62 = a;2 — 1 as functions of r are, in this regime, to first order, 

3 

= ei - ei + -ei 

= 62-62-262. (141) 

This corresponds to a focus in the (cuijCUi) plane, and hence with cui^h suffi- 
ciently small, 61 will have at least no zeros. □ 

Lemma VII. 5 If {coi,h,^^2,h) leads to a regular black hole solution with cui 
having fix nodes and lo^ having n^ nodes, then all {uJi^hi ^2,h) sufficiently close 
to {uJi^h, '^2,h) lead to regular or singular solutions with uui having at least fii 
zeros and L02 having at least n-i zeros. 

Proof 

Since the solutions are continuous in the starting parameters and r, for 
a;2,/i) sufficiently close to (cDi,^, a)2,fe), the gauge function lo\ will have 
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rii zeros and 002 will have n2 zeros for r < ri for some ri and uji{ri) will be 
close to their asymptotic values. □ 

From the proof of proposition |V11.1| , it is expected that the {LJi^h,^2,h) 
plane will be partitioned into open sets containing singular solutions with 
(ni,n2) nodes by lines of regular solutions. These lines of regular solutions 
will occur as the horizon radius, is varied. Numerical investigationsi'i'0 
indicate that for fixed Vh regular solutions exist at discrete points in the 
{uji^h,^2,h) plane, the positions of these discrete points varying as is var- 



ied. Lemma [VII. 5| shows that lines corresponding to regular solutions with 
different numbers of nodes must be disjoint, and furthermore that the region 
corresponding to singular solutions with (ni,?7.2) nodes will be bounded by 
lines of regular solutions having (ni,n2) nodes, where rii < ni and n2 < n2. 
Lemma |V11.5| is somewhat weaker than the corresponding result for su(2),l 



since the coupling between the cj's means that we cannot bound from above 
the number of zeros of the gauge field functions by analytic arguments. Be- 
low we shall address this question further by numerical investigations. This 



leads directly to our existence theorem [VII. 3| being slightly weaker than for 
= 2. 



Lemma VII. 6 Suppose (0,C(j2,h) corresponds to a charged regular solution, 
where UJ2 has n2 nodes. Given tiq, let Dq he the largest neighbourhood of 
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(0, uj2,h) for which {uji^h, ^2,h) correspond to regular or singular solutions with 
ui > uq and n2 > n2- Suppose, in addition, that there is (a)!,/!, a;2,/i) ^ Dq 
corresponding to a singular solution with ui = uq and n2 = n2- Then there 
is {uJ^i^,uJ2h) G Do corresponding to a regular solution with ni = uq and 
n2 = n2. 

Proof 

Let ho > no- Then there is a neighbourhood Do C Do, Do 7^ -Do for which 
all solutions are regular or singular solutions with ni > ho and n2 > n2- In 
Do there must be at least one singular solution having ni > ho and n2 > n2, 
by lemma [VII.5| , corresponding to starting values {oJi^h, ijJ2,h) ■ Consider a 



curve joining {LJi^h,^2,h) and (cDi,/i, ii'2,/1) and lying in Do- Then there must 
be at least one regular solution along this curve. Let {uj\,U2h) be the 
regular solution closest to {ui^h,'^2,h), having ni = nf and n2 = n^. Since 
{uj^f^,uj2 h) G Do, it follows that nf > no and > n2. Also, from lemma 
VIL5| , nf < no and < h2, since sufficiently close to i^2^h) there are 



singular solutions having Ui = no and ri2 = ^2. Therefore nf = no and 

nf = 77,2. □ 



Lemma VII. 7 Given n,2 and no > n2, then there exist regular and singular 
solutions with n2 = n2 and n\ > uq. 

Proof 
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Consider the point (0, uj2,h) which corresponds to the charged regular solution 
with ri2 = n2- Then all {uji^h,^^2,h) sufficiently close to (0,ti)2,h) are regular 
or singular solutions with Ui > uq and n2 > Let all such (wi /j, a;2,/i) 
form a neighbourhood Dq of (0,ci}2,/i)- From,! there exist (0,a)2,/i) £ -Do cor- 
responding to singular solutions having = and n2 = n2- By continuity, 
all ( i^2,/i) sufficiently close to (0,a;2,/i) correspond to singular solutions 



having ni > hq and 77-2 = n2. Hence, by lemma VII.6, there are also in this 



neighbourhood regular black hole solutions having rii > uq and 77.2 = ?t,2. □ 
Proof of Theorem |V11.3| 

Fix 77,2- Then we have regular solutions with node structure (n2,n2)- Now 
let hq =77,2 + 1, then from lemma |V11.7| there are regular solutions having 
7i2 = 77,2 and 77i > 77o. Let Uq be the smallest such ni. Now set no = Uq + 1 
and repeat the process. □ 

In order to guarantee the existence of black hole solutions having node 
structure (77,1,712) for every pair of integers (771,722), we would require the 
following lemma in addition to lemma |V11.7| . 

Lemma VII. 8 Suppose (cji /j, ci}2,h) corresponds to a regular black hole solu- 
tion in which Ui has fii nodes and UJ2 has fi2 nodes. Then, for {uji^h-,^2,h) 
sufficiently close to {ui^h, ^^2,h) , solutions for which UJ2 still has 77,2 nodes are 
such that uii has either 77,1 or 77,1 + 1 nodes. 

46 



The argument we now give for lemma |V11.8| is not an analytic proof 
because we require a numerical analysis. 

Proof 

By continuity, we know that for {u!i^h,^2,h) sufficiently close to {LJi h,LJ2,h), 
the gauge function ui has ni zeros and uj2 has n2 zeros for r < ri for some ri 
and uJi{ri) will be close to their asymptotic values. Since we are considering 
only solutions for which UJ2 has n2 nodes, for r > ri, then, 002 will be of one 
sign. Therefore we may consider the new dependent variable 

ij='^ (142) 

i02 

which will have the same number of zeros as c<Ji(r) for r > ri. Then the 
equation satisfied by ip is 

uj' 3 

Pfi^" + (2m - 2r%) ^' + 2rV— ^' + ^z^jj^ (l - = 0. (143) 

On some interval, ri < r < r2, the geometry will be very nearly fiat, and it 
will remain fiat until one of uJi, UJ2 or their derivatives approach O(y^). Next 
consider the fiat space equations for ip and a;2, which have the autonomous 
form (where r = logr): 

= U}2 - UJ2 + - + ^ii^ul^ UJ2 

= + —ijj + -ujip (1 - tp^) . (144) 

UJ2 2 
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Without loss of generality, we may assume that we are close to the critical 
point where tp = 1 and uj2 = V^- Using the substitution 

^ = l + ei, u;2 = ^2(1 + 62), (145) 

the linearized equations close to this critical point are: 

= e'l — ei — 6ei 

= es - £2 + 2ei - 2e2. (146) 

The analysis of section has already shown that we have a saddle point here; 
this notation is convenient merely because the perturbation in ip decouples 
from that of UJ2 in the linear approximation. From the form of the equation 
( |143| ), it is clear that will be monotonic increasing for > 1. In other 
words, for ei > initially, ei > always and ooi will have no additional zeros 
for r > ri. Hence we need only consider the case where ei < initially. 
Unfortunately the non-linear perturbation equations cannot be integrated 
analytically and so we present a numerical argument. For initial 61,62 suffi- 
ciently small, the values of ei and €2 will remain close to those on the unstable 
manifold for all t < ti for some ti, where Ti can be taken to be as large as 
we like by making the initial perturbations sufficiently small. Numerical in- 
tegration of the non-linear equations with the initial point on the unstable 
manifold and ei < 0, shows that €2 increases monotonically and ei decreases 
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monotonically with just one zero, and cuts through ei = —2 (corresponding 
to if) = —1), from whence ip must be monotonically decreasing (see figures 4 
and 5). Hence in this case Ui has ni + 1 zeros. □ 



Lemma [V 11 . 8| allows us to prove that every combination of integers (ni, 77-2) 



must correspond to a black hole solution. In the proof of theorem |V11.3 



we begin with regular solutions with node structure (n2,n2) and then from 



lemma [VII. 8| there are either regular or singular solutions with node structure 



n2 + 1, ^2). Using lemma [VII. 6| , there are then regular solutions having this 



node structure and the proof of theorem VII. 3 follows as before. 



So far we have proved the existence of infinitely many su(3) black holes 
only for rh > 2. The remainder of this section will be spent proving the result 
for < 2. In this case, the points at which the map / is not continuous 
(which must still exist by the argument used in proving proposition [VII. 1[) 
do not necessarily correspond to regular black hole solutions: they could also 
be oscillating solutions. A couple of lemmas concerning oscillating solutions 
are required before the proofs of proposition [VII. 1| and theorem [VII. 3[ can be 
extended. 

Lemma VII. 9 If {pi^h-,^2,h) leads to an oscillating solution with uji{t) 
anduj2{T) 1 asT ^ 00, anduj2 hasn2 zeros, then all {uJi^hi^2,h) sufficiently 
close to {oui^h, <^2,h) lead to one of the following types of solution: 
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1. an oscillating solution in which uji{t) and uj2{t) ^ 1 as t ^ oo, 
with UJ2 having at least n2 zeros; 

2. given uq, either a regular or a singular solution with Ui > uq and 
n2 > n2. 

3. an Soo solution. 
Proof 

Since in the original solution ui has infinitely many zeros and the solutions 
are analytic in r and the starting values, there is some Ti such that all 
solutions with (LJi^h,^2,h) sufficiently close to ci)2,/i) must have at least 
uq zeros of Ui and n2 zeros of 002 for r < ri, and all variables will be close to 
their asymptotic values. □ 

Lemma VII. 10 There exists {uji^h,^2,h) G D corresponding to an oscillating 
solution. 

Proof 

From,@ we know that there are points on the line uji^h = ^2,h corresponding to 
Soo solutions. By continuity, there is a neighbourhood Dq of the line uji^h = 
^2,h corresponding to Soo solutions. However, lemma |V11.5| still applies in this 



case, so that the regular solutions on the line uji^h = ^2,h will have singular 
solutions in D sufficiently close to them for which r(ro) > 2. Therefore, 
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by lemma |V11.9|, there must be points in D which correspond to oscillating 



solutions. □ 



The proofs of proposition |V11.1| and theorem |V11.3| are now exactly the 



same as for rh> 2. Since oscillating solutions have an infinite number of zeros 
for at least one of the gauge field functions, and using lemma |V11.9|, regular 



solutions are still needed to separate singular solutions having different node 
structures. 

VIII. su(iV) black holes 

The detailed discussion of the previous section now enables us to proceed 



quite swiftly to the analogues of proposition [VII. 1| and theorem [VII. 3| for gen- 
eral N . The method of proof will be by induction, which was the basic idea 
used in the previous section, where we proved existence for su(3) exploiting 
known results about su(2). 

The first step is to illustrate how solutions which are from su(n), where 
n < N may be embedded in the su{N) framework. Firstly, exactly as in 
the AT = 3 case, we may embed neutral .u(2) solutions.0 We set c.,(r) = 
- j)uj{r) and introduce a scaled variable R = X^r, where 

A^= QiV(iV + l)(iV-l)y . (147) 
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Then the field equations are exactly the same as the su(2) ones, and the 
existence of solutions follows directly .0 Secondly, charged, effectively su{N— 
1) solutions can be generated by setting ui = (or ujj^^i = 0).i Although 
the equation for m' in this situation is not the same as for neutral su{N — 1) 
black holes, due to the charge, the regular equations discussed in section || 
are unchanged, so the proof of existence of neutral 3u{N — 1) black holes 
extends naturally to this case. 

For > 3, there are additional embeddings which arise from setting at 
least one of uj2, ■ ■ ■ , ujn-2 = 0, so that the gauge field equations decouple into 
two or more sets of coupled components, the sets being coupled to each other 
only through the metric. Solutions of this form have been found numerically 
in.i We illustrate how the existence of black holes of this type may be proved 
by considering the simplest case, which arises when = 4. In this case, we 
have three non-zero gauge field functions, uji,uj2,ujs. If we set UJ2 = 0, then 
the field equations take the form 




(2m - 2r^pe) uj[ - [l - uf) Ui 



(148) 




(2m - 2r^pe) - (l - uj) UJ3 



(149) 



m 



(150) 



2G 



r 



(151) 
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where 



G = co[' + u'l = ^ (K - 1]% H - 1]' + 8) . (152) 

These equations look very much hke two uncoupled su(2) degrees of freedom, 
however, the two cu's are (albeit weakly) coupled. The fact that we have two 



gauge field functions here means that we can use the methods of section [VII 
to prove the existence of regular black hole solutions. In fact, all the results 
of that section carry directly over to this situation on replacing uj2 there by tus 



here. There is one exception. The equations ( |148| , |149| ) are slightly different 



from those in section |VI1| and this enables us to strengthen lemma |V II. 5| . 

Lemma VIII. 1 // (cji jj, 0)3 jj) leads to a regular black hole solution with ui 
having ni nodes andco^ having nodes, then all {(jji^h,^z,h) sufficiently close 
to {ijJi^hi^z,h) l^o-d to regular or singular solutions with uji having either fii 
or ni + 1 zeros and uj^ having either or + 1 zeros. 

Proof 

Since the solutions are continuous in the starting parameters and r, for 
(tui /i, ciJs /j) sufficiently close to (c^i /j, 0)3^/1), there is some ri such that uJi will 
have fii zeros and ^^3 will have zeros for r < ri. There will also be an 
r2 > ri such that all field variables are close to their asymptotic values for 
r G [ri, r2]. The fiat space field equations (see section |Vp decouple completely 
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since . 0, and so the phase plane analysis ol is valid in this case. The 
phase plane analysis enables us to be more precise about the number of zeros 
of the perturbed gauge field functions, unlike the more general case where 
the phase space had more than two coupled dimensions. Therefore the con- 
clusions o^ can be applied directly here and each gauge field function can 
have at most one zero for r > r2- □ 



With this more powerful lemma, the analysis of section |V11| reaches the 



conclusion that there are regular black hole solutions of the required form 
for each rh and node structure (ni, ns). 

For general N, similar embeddings of the form of two or more su(n) 
{n < N) type solutions separated by at least one loj = are possible, and the 
existence proof follows the lines outlined above for su(4), using the existence 
of solutions for the various su{n) . The allowed node structures will be exactly 
the same as those for the constituent su{n) components. We refer the reader 
to@ for further details of the construction of solutions of this form. 

The remainder of this section will be spent outlining briefly the proof of 
the following theorem. By 'types' in the statement of the theorem we mean 
the node structures of the gauge field functions. 

Theorem VIII. 2 There exist infinitely many types of regular black hole so- 
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lutions of the su{N) Einstein- Yang-Mills equations including genuine su{N) 
solutions and embedded solutions. 



All that remains to prove this theorem is to show the existence of genuine 
su{N) solutions in which the node structures of the gauge field functions are 
not all the same. We assume for an inductive hypothesis that the theorem 
holds for all n < N. The situation is slightly complicated by the fact that 
our parameter space consisting of the Wj^hS is now — 1-dimensional. We 
consider the uji^h = hyperplane, on which we know the solution space by 



the inductive hypothesis. The map / (see section [VI 1|) now takes M^-i to 



^ X [0, 1]. Again, we need only consider the action of / on those parts of 
iN-i which no ujj^h vanishes. We restrict attention to one of the disjoint 



spaces so generated without loss of generality. Lemma [VII. 2| and the proof 
of the su{N) equivalent of proposition [VII. 1[ now carry over to prove the 
existence of genuine su(A^) black holes. Lemmas [VII. 4[ , |VII.6[ and |VII.7[ now 
hold for a point (0, LJ2,h, ■ ■ ■ , ^n-ia) corresponding to a charged solution, and 
lemmas [VII. 5| and [VII. 9[ carry straight over to a point {cui^h, ^2,h, ■ ■ ■ , ^N-i,h) 
which leads to a regular black hole solution or oscillating solution respectively. 
These are the ingredients necessary for the simple argument which leads to 
the proof of theorem [VII. 3[ and hence theorem [VIII. 2[ for both r\ > \N[N + 
1)(A^ — 1) (in which case no oscillating solutions exist), and < \N{N + 
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1)(A^— 1). Note that, as in the su(3) case, we are not able to prove analytically 
that every sequence of integers (ni, 77-2, . . . , '"-Af-i) corresponds to a black hole 
solution whose gauge functions have this node structure, although it might 
reasonably be expected that this is indeed the case. 

IX. Results and conclusions 

In this paper we have proved the existence of a vast number of hairy black 
holes in su{N) Einstein- Yang-Mills theories. The solutions are described by 
N-1 parameters, corresponding to the number of nodes of the gauge field 
functions. The result of0 tells us that each of these solutions will have 
a topological instability, similar to the flat-space "sphaleron", as was the 
case for su(2) black holes. This instability does not necessarily diminish the 
physical importance of these objects, as they may have an important role in 
processes such as cosmological particle creation.lli 

We now briefly mention another important reason for studying su{N) 
black holes, namely the behaviour of the solutions as — 00. Black holes 
in su(oo) Einstein- Yang-Mills theory would possess an infinite amount of hair, 
requiring an infinite number of parameters to describe their geometry. This 
might be analogous to the "W-hair" found in non-critical string theory ,@ and 
have drastic consequences for the Hawking radiation, information loss and 
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quantum decoherence processes associated with black holes. The existence 
of infinite amounts of hair might also render such objects stable. It is already 
known that the su(cxd) Lie algebra is simply that of the diffeomorphisms of 
the sphere. This means that the limit as ^ oo cannot be taken smoothly, 
in particular the results of the present paper are only valid when is finite. 
We hope to return to this matter in a subsequent publication. 
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Figure 1: The (<^i,/i, <^2,/i) plane, with a schematic representation of the hnes 
representing regular black hole solutions, which are also denoted by crosses. 
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Rm — 1 



Rm — 



A^i = 



A^i = 0.5 



Figure 2: Schematic representation of the set B x [0, 1], corresponding to the 
space of parameters of the su(3) solutions in section VII. The horizontal axis 
represents A^i, the number of nodes of the gauge function cui, and the vertical 
axis Rm, which characterizes the maximum value of the radial co-ordinate r 
for each solution. 
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^ (g) (2)(a88) ® 



Figure 3: The projection of the parameter space of solutions on to the set 
B X B. The axes are A^i and N2, the number of nodes of the respective gauge 
functions uji and UJ2. The crosses indicate schematically the points in this 
set. 
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-0.5-- 



e[l] (tau) 




-1.5 



Figure 4: The su(3) perturbation ei{r) along the unstable manifold. The 
value T — corresponds to the critical point. Note that ip — 1 + ei mono- 
tonically decreases through and —1, from which point it will continue to 
monotonically decrease. 
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e [2] (tau) 



1-- 



0.5-- 



tau 



Figure 5: The su(3) perturbation 62 along the unstable manifold. Note that 
this is a monotonically increasing function, which is in accordance with the 
gauge function lo^ having no further zeros. 
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